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ABSTRACT 

We present a multi- wavelength photometric study of ~ 15, 000 resolved stars in the nearby spiral 
galaxy M83 (NGC5236, L> = 4.61Mpc) based on Hubble Space Telescope Wide Field Camera 3 ob- 
servations using four filters: F336W, F438W, F555W, and F814W. We select 50 regions (an average size 
of 260 pc by 280 pc) in the spiral arm and inter-arm areas of M83, and determine the age distribution 
of the luminous stellar populations in each region. This is accomplished by correcting for extinction 
towards each individual star by comparing its colors with predictions from stellar isochrones. We com- 
pare the resulting luminosity weighted mean ages of the luminous stars in the 50 regions with those 
determined from several independent methods, including the number ratio of red-to-blue supergiants, 
morphological appearance of the regions, surface brightness fluctuations, and the ages of clusters in 
the regions. We find reasonably good agreement between these methods. We also find that young 
stars are much more likely to be found in concentrated aggregates along spiral arms, while older stars 
are more dispersed. These results are consistent with the scenario that star formation is associated 
with the spiral arms, and stars form primarily in star clusters and then disperse on short timescales 
to form the field population. The locations of Wolf-Rayet stars are found to correlate with the po- 
sitions of many of the youngest regions, providing additional support for our ability to accurately 
estimate ages. We address the effects of spatial resolution on the measured colors, magnitudes, and 
age estimates. While individual stars can occasionally show measurable differences in the colors and 
magnitudes, the age estimates for entire regions are only slightly affected. 
Subject headings: galaxies: individual (M83, NGC 5236) — galaxies: stellar content 
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1. INTRODUCTION 

Understanding the properties of stars and the history 
of star formation in galaxies remains one of the most 
fundamental subjects in astrophysics. The Hubble Space 
Telescope (HST) provides an important tool for this en- 
deavor, since it enables the detailed study of stars and 
star clusters, not only in the Milky Way and its nearest 
neighbors, but in galaxies well beyond the Local Group. 
The Wide Field Camera 3 (WFC3) provides a particu- 
larly valuable new panchromatic imaging capability, with 
spectral coverage from the near-UV to the near-IR. This 
capability is especially useful for studying the stellar pop- 
ulations in nearby galaxies where individual stars are re- 
solved. 

A good example of a multi- wavelength survey of nearby 
galaxies which employs resolved stars and star clusters 
is the ACS Nearby Gal axy Treasury (ANGST) program 
( |Dalcanton et al.][2QQ9 ). This study includes ~65 galax- 
ies out to ~3.5 Mpc, and provides uniform multi-color 
(BVI) catalogs of tens of thousands of individual stars 
in each galaxy. Although the ANGST program provides 
an excellent survey for a wide range of studies, it does 
not provide observations in the U -band since the observa- 
tions were obtained before WFC3 was installed on HST. 
The L^-band is particularly useful for age-dating popula- 
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tions of young stars, which is the focus of the current 
paper. PHAT (Panchromatic Hubble Andromeda Trea- 
sury) is a related study that does take advantage of the 
new /7-band capability on WFC3. While nowhere as ex- 
tensive as the Multiple Cycle Treasury Program PHAT 
survey, the current study of M83 is complementary in the 
sense that it provides similar observations of the resolved 
stellar component of a nearby spiral galaxy, and uses 
quite different analysis techniques, as will be discussed in 
Sections 2 and 3. Future observations of M83 (ID:12513, 
PI: William Blair) will expand the WFC3 dataset avail- 
able for M83 from 2 to 7 fields. 

M83 (NGC5236), also known as the "Southern Pin- 
wheel" galaxy, is a slightly barred spiral galaxy with a 
starbursting nucleus located at a distance of 4.6 1 Mpc, 
i.e., (m - M) = 28.32 ± 0.13 ( |Saha et al.pOO^ . The 
Ha emission can be used to pinpoint regions of recently 
formed massive stars along the spiral arms, while red su- 
pergiants can be found throughout the galaxy. The color- 
magnitude diagram (CMD) of resolved stars is a power- 
ful diagnostic tool for understanding the stellar evolution 
and history of star formation of galaxies in detail. By 
comparing stellar evolution models to observed CMDs, 
we are able to determine the ages of the stellar popula- 
tions in galaxies. However, in active star-forming regions, 
the stars are often partially obscured by dust. Applying 
a single extinction correction value for the whole galaxy 
often results in over- or under-estimates of the ages of in- 
dividual stars. If the CMD is the only tool used to deter- 
mine the ages, the spatial variations of dust extinction in 
galaxies and their effects on the determined ages are not 
readily apparent. The additional information available 
from the color-color diagram can remedy this problem. 
By using techniques developed in this paper, we can con- 
strain the variations of dust extinction across M83, and 
make corrections for individual stars. We also focus on 
spatial variations in the stellar properties throughout the 
galaxy. 

The HST WFC3 observations of M83 were performed 
in August 2009 as part of the WFC3 Science Oversight 
Committee (SOC) Early Release Science (ERS) program 
(ID: 11360, PI: Robert O'Connell). The central region 
(3.2kpcx3.2kpc) of M83 was observed in August, 2009. 
A second adjacent field to the NNW was observed in 
March 2010, and will be included in future publications. 
Details of the W FC3 ERS data calibrat ion a nd process- 
ing are given by Chandar et al. | ( |2010[ ) and |Windhorst] 
et gllpQTT ). 

In this paper we present UBVI photometry of re- 
solved stars in M83 and the resulting color-color and 
color-magnitude diagrams. We use (F336W — F438W) 
vs. (F555W — F814W) color-color diagrams to con- 
strain the variation of dust extinction along the line- 
of-sight of each individual star, and correct FS14W vs. 
(F555W-F8UW) CMDs for the extinction of individual 
stars to determine their ages. The high sensitivity and 
the superb resolving power in the WFC3 F336W-band 
plays a key role in allowing us to develop our extinction 
correction techniques, and demonstrate the performance 
of the newly installed WFC3. This results in improved 
stellar age estimates, and a better understanding of the 
recent star- formation history of M83. 

Our investigations are focused on the followings: 
(1) Do we see spatial variations of stellar ages in M83? If 



so, can we use these variations to learn more about 
the evolution of the galaxy and what triggers star 
formation? 

(2) How well do various age estimates correlate (i.e., re- 
solved stars, integrated light, clusters, number ra- 
tio of red-to-blue stars, Ha morphology, stellar sur- 
face brightness fluctuations, presence of Wolf-Rayet 
stars)? 

This paper is organized as follows. In §2, we describe 
the observations, photometric analysis, and the extinc- 
tion correction method. CMDs and color-color diagrams 
of the 50 selected regions are then used to provide stel- 
lar age estimates, as described in §3. Comparisons with 
other age estimates (i.e., integrated light and star clus- 
ters) are also made. In §4, we compare the stellar ages 
to a variety of other parameters that correlate with age, 
including red-to-blue star ratios, Ha morphology, and 
surface brightness fluctuations. Section 5 includes a dis- 
cussion of how these comparisons can be used as diagnos- 
tics, with special attention paid to the question of what 
spatial variations can tell us about the star formation 
history in M83, and on the question of how star clus- 
ters might dissolve and populate the field. A summary 
of our primary results is provided in §6. In Appendix 
A, we investigate the locations of sources of Wolf-Rayet 
stars in our M83 field and discuss the correlation with 
the positions of young regions. In Appendix B, we de- 
scribe how the spatial resolution affects our measured 
colors and magnitudes, and our age estimates of stars in 
our M83 field. 

2. HST WFC3 OBSERVATIONS 

2.1. The WFC3/UVIS Data 

The data were obtained in seven broad- and seven 
narrow-band filters in the WFC3 UVIS and IR chan- 
nels. Details are describe d in Dopita et al. (2010) 
and |Chandar et al 
use tour broad- ban* 



7] ( |20l0l ). 

a images. 



In the current paper we 
The filters and exposure 
times are F336W (1890 sec), F438W (1880 sec), F555W 
(1203 sec), and F814W (1203 sec). Although we do 
not make transformations to the Johnson-Kron-Cousins 
UBVI system, we will refer to these filters as J7, V, 
and / for convenience. Three exposures were taken in 
each filter at different dithered positions to remove cos- 
mic rays, the intrachip gap, and to partially compensate 
for the undersampled point-spread function (PSF) of the 
WFC3/UVIS channel. The raw data were processed us- 
ing the MULTIDRIZZLE software flKoekemoer et al. [2002 ) 
with an effective pixel size of / /0396, which corres ponds 
to 0.885 pc per pixel at the distance of 4.61 Mpc (Saha 
|et al. |2006[ ) . The resulting multidrizzled images are com- 
bined, cosmic-ray removed, and distortion-corrected. A 
color composite of the F336W, F555W, F814W (broad 
bands) and F502N, F657N (narrow bands) images is 
shown in Figure [I] It covers the nuclear region of M83, 
part of its eastern spiral arm, and the inter-arm region. 
Stars in the cores of compact star clusters are not re- 
solved at this resolution, but the brighter stars in the 
outskirts of the star clusters and in the field are gener- 
ally resolved into individual stars. Details of the effects of 
spatial resolution are discussed in Appendix B. Since we 
are interested in spatial variations of stellar ages in M83, 
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we selected 50 regions in the spiral arm and inter-arm ar- 
eas. Boxes outlined in blue (very young - age < 10 Myr, 
see §3), yellow (young - 10 < age < 20 Myr), and red 
(intermediate- aged - age > 20 Myr) in Figure [l] show the 
50 regions selected for detailed study in this paper. 



2.2. Photometry and Artificial Star Tests 

Photometric analysis of the WFC3 M83 data was per- 
formed o n the U, B, V, and I images using the DoPHOT 
package (Sche chter et al. |1993 [) with modifications made 
by A. SaEal Additional routines to derive parameters 
for the kurtosis of the analytic PSFs used by DoPHOT, 
and for post-processing DoPHOT output to obtain cali- 
brated aperture-corrected photometry, were performed 
using customized IDL code written by Saha. A more de- 
tailed descr iption of these procedures can be found in 
|Saha et al. (2010). The kurtosis terms were derived by 
fitting the functional form of the DoPHOT analytic PSF 



tting the 

to 20 relatively bright and isolated stars. The remain- 
ing shape parameters (a x , a y , and a xy ) are dynamically 
optimized within DoPHOT. 

The difference between the shape parameters for an in- 
dividual object and those for a typical star was calculated 
and used to classify the object as a star, a galaxy, or a 
double-star. The shape parameters were determined sep- 
arately for each image. The numbers of objects classified 
as stars are approximately 20,000 in U and £?, 17,000 
in V, and 15,000 in /band images, respectively. The 
calculated Full Width at Half Maximum (FWHMs) are 
U = 0'.'096 (2.42 pixel), B = 0'/093 (2.34 pixel), V = Of! 091 
(2.29 pixel), and / = 0'/113 (2.85 pixel). All magnitudes 
are in the WFC3-UVIS VEGA MAG magnitude sys tem, 
calculated using equation (4) in |Sirianni et al. | ( |2005| ) and 
the latest zeropoint magnitudes! U = 23.46, B = 24.98, 
V = 25.81, and 7 = 24.67 mag, provided by STScI at the 
HST/WFC3 website 22 . 

A multi- wavelength UBVI catalog was then con- 
structed by matching the individual objects in each cat- 
alog between these filters. Interestingly, we found that 
the extreme image crowding in these images, combined 
with the excellent panchromatic sensitivity and unprece- 
dented spatial resolution provided by the newly installed 
WFC3, introduced a new challenge during this matching 
step, with large numbers of "pseudo" matches occurring 
if too-large a match radius was used. We discuss the 
effects of spatial resolution in Appendix B. 

Figure [2] demonstrates this problem by showing the 
F555W and color composite image cut-outs near Regions 
# 26 and # 27. In the middle panel of Figure [2| stars 
circled in red are objects with (V—I) colors that are very 
red ((V — I) > 1.2 mag), but have (U — B) values that 
are fairly blue (—0.4 < (U — B) < —1.7 mag). These 
can be seen as the spray of points in the upper right 
parts of the top right panels of Figures [3] and [4j which 
will be discussed below. As shown by the postage stamp 
images in Figure [2| about half of these are likely to be 
blue stars with very high values of reddening (i.e., the 
bottom panel and the small postage stamps on the right 
show that they are in dusty regions) , while the other half 
are very close superpositions of at least one red and one 
blue star, which we will call "pseudo" matches (i.e., the 



panels on the left show a strong blue and red gradient 
across many of the objects). 

We found that using a very small matching radius of 
0.5 pixels was needed to minimize "pseudo" matches in 
our catalog. Earlier attempts using a matching radius 
of 3 pixels resulted in ~10 % pseudo matches. Accu- 
rate matching also requires precise geometric distortion 
corrections in all filters. However, even when these con- 
ditions are met, there is a finite chance that two different 
stars will fall within the same aperture. This is most eas- 
ily seen by blinking the U- and /-band images. While the 
single stars stay in the same position, the pseudo match 
stars move slightly, showing they are two different stars 
with different colors. 

The final stellar catalog contains ~15,000 objects. We 
find that our procedure results in only a few percent of 
pseudo matches, as will be discussed in §2.4 (i.e., the 
upper right panel of Figure [4]). 

To measure the photometric completeness, we per- 
formed artificial star tests with the same detection and 
photometry procedure as applied to the actual stars. We 
inserted 1000 artificial stars with Gaussian PSFs and 
FWHMs appropriate for these stars at random positions 
into all four images. The magnitude of the inserted stars 
varied from 22 to 28 mag in steps of 0.25 mag. The 
50% completeness levels at 5 a detection thresholds for 
a typical region are reached at recovery magnitudes of 
approximately 24.3, 25.1, 25.3, and 24.8 mag for the U, 
B, V, and / images, respectively. These limits can vary 
by about a magnitude, depending on the brightness of 
the background in a given region. The exception is Re- 
gion #48, which includes the nucleus. The completeness 
thresholds are two or more magnitudes brighter in this 
region due to the very high background. This region has 
therefore been excluded from the rest of the discussion, 
but is included in Table Q] for reference. 

2.3. CMDs and Color-color Diagrams in M83 

The top panels of Figure [3] show the / vs. (V — I) CMD 
(left) and the (U — B) vs. (V — /) color-color diagram 
(right) of resolved stars extracted from the cross-matched 
catalog. The CMD shows the presence of young main- 
sequence (MS) stars, transition or He-burning "blue 
loop" stars, red giant stars (hydrogen shell-burning and 
hydrogen+helium shell-burning stars), and a relatively 
small number of low-mass red-giant branch (RGB) stars. 
This is because t he tip of the low-ma ss RGB stars is 



at /= 24.7 mag ( Karachentsev et al. 2007[ ), which is 



roughly the same as our completeness threshold. The 

100 Myr. The 
ranging in age 



ages of these stars range from ^ 1 
Padova isochrones (Marigo et al. 
from log r (age/yr) 



to 



2008D 

6.55 (3.51viyr) to log r (age/yr) 



8.0 (100 Myr) are included in Figure [3] for reference. 
We note that the younger isochrones (e.g., 1 Myr) fall 
nearly on top of the 3.5 Myr isochrone, and hence are 
not included in the diagram. 

The arrows in both panels show Ay = 1 reddening vec- 
tors for M83 using t he Ry = 3.1 extinction curve of 
Cardelli et al. | ([1989 ) . Correct ions have been made for 
foreground Galactic extinction (Schlegel et al. 1998]) us- 
ing A F336W 0.361, A F438W = U.29U, A F555 w = 0.229, 
and A F8 i4w = 0.133 mag, respectively. 23 In Figure [3) 



22 http://www.stsci.edu/hst/wfc3/ 



23 NASA/IPAC Extragalactic Database (NED): 
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we only plot stars with photometric errors less than 0.25 
mag in the filters used for each diagram. The numbers of 
objects plotted in the CMD and the color-color diagram 
are about 12,000 and 8,500 stars, respectively. Since ap- 
proximately 30% of the stars are not detected in one or 
more of the ?7, B, or / filters, due to either reddening or 
intrinsically red or blue colors, fewer stars are plotted in 
the color-color diagram. 

We adopt a value of 1.5 times solar metallicity 
(Z=0.03), the highest metallicity isochrones available 
from the Padova dat abase, for M83 based on Bresolin 
fc Kennicutt | ( |2002[ ). A test using solar metallicity 
(Z=0.019j isochrones showed that the ages for a typical 
region would be ~ 2 Myr older if the low er metallicity 
i s ado pted. This agrees with the study by Lar sen et al. | 
(2011), who found that the simulated CMDs of young 
massive clusters in M83 with solar and super-solar metal- 
licity isochrones would not look much different. The- 
oretical isochrones calculated for the WFC3 filters and 



Z=0.Q3 (1. 5 Z ) from the Padova database 24 ( |Marigo 
et al. |2008| ) are overlaid on both panels in Figure [3] The 
dashed line in cyan shows the 50% completeness level i 



cyan snows tiie 5U7o completeness level m 
/ and (V—/), using the completeness threshold numbers 
from the artificial star test in §2.2. 

The bottom panels of Figures [3] and [4] show the his- 
tograms of ages and masses for the same stars plotted in 
the CMD and color-color diagram. The ages and masses 
of stars were estimated by finding the closest match be- 
tween the / and (V — /) values for each star, using a 
fine mesh of the stellar isochrones plotted in the CMD 
and color-color diagram in Figures [3] and [4j More de- 
tails about the stellar age-dating from the CMDs and 
the stellar isochrones are given in §3.1. 

2.4. Extinction Corrections 

Since M83 has an intricate structure of dust lanes asso- 
ciated with active star forming regions in the spiral arms 
and thin layers of dust in the inter-arm area in Figure [TJ 
we cannot use a single value of internal extinction and 
apply it to correct for the extinction of all stars in a given 
region. In the color-color diagram in Figure [3J we notice 
that while the bluest stars match the model isochrones 
quite well after we correct for Galactic foreground ex- 
tinction, the majority of stars are found redward of the 
models, implying a typical extinction of Ay ~ 0.5 mag 
(for the densest part of the data "swarm"). As described 
below, we can determine the reddening of these stars if we 
assume that they are intrinsically blue stars that belong 
on the Padova models (i.e., with (V — I) c± —0.3 mag). 
A visual inspection supports this interpretation. Stars 
with observed values of (V — I) ~ —0.3 mag are found 
in areas with no obvious dust features surrounding them 
while stars with (V — I) > 0.5 mag are near dust fila- 
ments. 

Our basic approach to correct the colors (and lumi- 
nosities) of individual stars for the effects of extinction is 
to backtrack the position along the reddening vector for 
each data point, until we hit the stellar isochrones in the 
color-color diagram. What makes this method work is 
the fact that all the isochrones are nearly on top of each 
other in the range (U — B) <0.0 mag. The way we apply 

http:/ /ned. ipac.caltech.edu/ 



this in practice is to match the observed and predicted 
values of the reddening free Q parameter defined by: 

Qubvi = {U-B)- ^yZf) X ( V - J ) (!) 



(Mihalas & Binney 1981|) . Using observed values of 
UBV1, we compute the observed Qubvi values by us- 
ing the standard slope of E(U — B)/E(V — I) = 0.58 
flMihalas fc Binney |[l98Tj | Whit more et al ||1999[) , which 
corresponds to Ky = 3.1. The predicted Q values are 
calculated using the 1.5 Z Padova isochrones. A compli- 
cation is that we cannot simply use the matched Q values 
to determine corrected values of (U—B) and (V—/), since 
they would, by definition, fall precisely on the isochrones 
(i.e., we would be using circular reasoning). Instead, we 
solve for the extinction in (U — I), which is partly inde- 
pendent and uses the longest wavelength baseline. We 
then calculate the extinction values in UBVI and the 
color excess values in (U — B) and (V — /). The averages 
of the corrected internal extinction and color excess val- 
ues by our star-by-star correction method are Apsmw = 
0.696, A F438W = 0.559, A FbbbW = 0.441, A F814W = 
0.256 mag, and E(U-B) = 0.137 and E(V-I) = 0.185 
mag. This is similar to the average internal extinction 
(ApbbbW = 0.425 mag) estimate d from the cluster SED 
fitting by |Bastian et al.' ( |2011[ ). However, our results 
differ from the average extinction (Ay = 0.671 mag) of 
45 clusters in the nucleus (~ 20" i n diameter) of M83 , 
determined from Ha/Hf3 ratios by Harris et al. | ( |2001[ ). 
This is consistent wit h the finding ol l arger extinction m 



CMD version 2.3; http://stev.oapd.inaf.it/cgi-bin/cmd 



the nuclear region bylChandar et al. I (120101). 

In the color-color diagram (top right panel) of Figure [3j 
we note that while most of the data points are consistent 
with the standard reddening vector, the stars located in 
the green triangle appear to follow a flatter reddening 
vector. Whether this is actually due to a different red- 
dening law in M83 (e.g., as has been suggested for heavily 
extincted regions such as # 4, # 12, and # 20), or is some 
sort of artifact (e.g. due to photometric uncertainties, 
unresolved star clusters, or the "pseudo" matching prob- 
lem discussed above) is difficult to determine from our 
present dataset. However, for our specific needs the pre- 
cise answer to this question is not critical. Our approach 
will be to correct the data points in the triangle back to 
a position near the top of the isochrones, as shown by 
the intersection of the red dotted line and the isochrones 
in Figure [4j This is equivalent to using a range in Ry 
between 3.1 and 5.7 (the value represented by the red 
dotted line). 

The upper panels of Figure [4] show the CMD and 
color-color diagram corrected for the internal extinction 
in M83. The arrows in both panels are the reddening 
vectors. We do not apply an extinction correction for 
stars with (V — I) > 1.2 mag, or above an extrapolation 
of the flatter extinction vector from the bluest possible 
isochrones (see the red dotted lines in the upper right 
panel of Figure [4|, since we believe that the colors of 
many of these sources are incorrect (i.e., roughly half 
are likely to be pseudo matches, as discussed in §2.2 and 
shown in Figure |2| . We note that only about 4 % of 
the stars fall in this part of the diagram. These stars 
are removed from the subsequent analysis. If we include 
these stars, our age estimates typically increase by about 
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lMyr. 

Also, stars below the blue dotted line in Figure [4] are 
not corrected for extinction, since these are main se- 
quence turnofT stars (i.e., blue loop or transition stars 
in general), and hence, do not have unique UBVI col- 
ors for the stars along a given reddening line, unlike the 
main sequence stars with (U— B)<0.0 mag. 

Based on the extinction corrected CMD in Figure [4j 
we note that: (1) the corrected data now shows a rel- 
atively narrow distribution of points along the left side 
of the CMD in agreement with ages of ~3 Myr, and 
(2) while the data now aligns with the model isochrones 
much better, there is still a fair amount of scatter (espe- 
cially for / < 24), primarily due to observational uncer- 
tainties from the fainter stars. 

The bottom panels of Figure [4] show the distribution 
of ages and masses of stars after we apply our extinction 
correction technique. Based on the automatic method 
described in §3.1, we note that the luminosity weighted 
mean ages determined from the corrected colors and lu- 
minosities are significantly younger (~14 Myr; Figure 
|4| than the mean ages determined from the uncorrected 
colors and luminosities (~27 Myr; Figure |3|, as we ex- 
pected. 

3. AGE-DATING POPULATIONS IN M83 

In this section we discuss several different methods of 
age-dating stellar populations in M83, using both indi- 
vidual stars and star clusters, and then intercompare the 
results. We do not necessarily expect field star and clus- 
ter ages to agree, since the dissolution of clusters may 
bias the observed cluster population towards younger 
ages than the surrounding field stars. 

3.1. Age-Dating the Resolved Stellar Population using 
Color Magnitude Diagrams 

As shown in Figure [TJ we selected 50 regions that cover 
spiral arm and inter-arm areas in order to study the re- 
cent star- formation history of M83. Figures [51 and [6] show 
cutouts of Regions # 12, # 7, # 29, and # 2?* These fig- 
ures include color-color diagrams of the stars in these re- 
gions, as well as CMDs that are uncorrected (upper) and 
corrected (lower) for extinction. These four regions high- 
light the range in age of the dominant stellar population, 
from very young to intermediate ages: i.e., Region #12 
with very strong Ha emission superposed on the cluster 
stars; Region # 7 with a large bubble of Ha emission 
surrounding the stars and clusters; Region # 29 with no 
Ha emission, but large numbers of bright red and blue 
stars; Region 2 with no Ha emission and fainter stars. 

A cursory glance at the CMDs of these four regions 
(Figures [5] and [6| show clear differences. The primary dif- 
ferences are, as predicted by the isochrones, that younger 
regions contain: (1) bluer main sequence stars; (2) larger 
numbers of upper main sequence stars; and (3) larger 
ratios of blue to red stars. Our regions do not, however, 
contain only stars of a single age. Even these four re- 
gions, which were chosen to include a single dominant 
stellar population, contain a mix of young and old stars. 

25 The color-magnitude diagrams and color-color diagrams 
of all 50 regions (similar to Figures 5 and 6) are avail- 
able on the WFC3 Early Release Science Archive website 
(http:/ /archive. stsci.edu/prepds/wfc3ers/ ). 



In this paper, we are primarily interested in the ages of 
the bright young stars in these regions, which dominate 
CMDs in luminosity-limited samples. 

In practice, this focus on the younger population is 
carried out by imposing a magnitude cutoff of Mj = 
—5.5 mag (i.e., / = 22.82 mag which corresponds to the 
age cutoff of ~60Myr in the CMD) for the stars used 
to estimate ages (i.e., the magenta dashed line in Figure 
[4|. This luminosity limit allows us to stay above the 
completeness limit for all but the very reddest stars (see 
Figure [3|, while also focusing on evolved stars and those 
on the upper main sequence, which are most sensitive 
to the age of a stellar population. Our primary goal is 
to obtain reliable "relative" estimates for the 50 regions 
(rather than absolute ages), and our magnitude limit is 
sufficiently deep to accomplish this goal. 

We compare the extinction corrected colors and mag- 
nitudes of stars with the Padova stellar isochrones in two 
ways to estimate the age of the dominant stellar popu- 
lation in each region. First, two of the authors (HK and 
BCW) independently estimated ages based on a visual 
comparison of the CMDs and isochrones. We focused 
on features such as the color of the main sequence stars, 
the number of stars in the upper main sequence, and 
the number ratio of blue to red stars. The left panel 
of Figure [7] shows that the independent, visually deter- 
mined age estimates are in good agreement with a 7.9 a 
(slope/uncertainty of the best linear fit) correlation and 
a slope within 1 a of the unity value. The average ages 
of these two independent (manual) estimates are listed 
in column Age C MDman in Table [I] 

Next, we estimated the ages automatically, finding 
the closest match in both age and mass for each star 
by comparing the extinction corrected / magnitude and 
(y — I) color with a fine mesh of stellar isochrones gen- 
erated from the Padova models (logr (age/yr) ranging 
from 6.05 to 8.35 in steps of 0.05). We then calculate 
a luminosity-weighted age from the distribution of in- 
dividual stellar ages in each region. The result of the 
automatic age estimate is listed in column AgecMDauto 
in Tabled] 

To assess the uncertainties in our age and mass deter- 
minations for each star, we ran a test for four different 
cases by adding and subtracting the photometric errors 
in the (V — I) color and /-band magnitude. This moves 
each object on the CMD to the right (+(V—I) error ), left 
(-(V -I) error ), down (+m /err , or ), and up (-m Ierror ). We 
then estimated the ages for these four cases in the same 
way as we did to get AgecMDauto- The results of this 
test are summarized in column AAgecMDauto in Table 
[I] As shown in the AAgecMDauto column in Table [I] 
the variation in the age estimates are typically only a 
few Myr (i.e., <10%). 

Figure [8] shows the age and mass histograms of indi- 
vidual stars in the four regions shown in Figure [5] and [6] 
Luminosity-weighted mean values of the stellar ages are 
shown in the upper right of the left panels. The age se- 
quence of these regions is as expected, from the youngest 
in Region # 12 to the oldest in Region # 2. 

One unexpected result was that while Region 2 was 
selected to be representative of an older field population, 
it also contains a fair number of bright young stars. This 
suggests that young individual stars may form in the field 
in low numbers. Another possibility is that many of these 
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are "runaway" stars, i.e., young stars that have been dy- 
namically ejected from their birth-sites within young star 
clusters due to interactions with other stars. IWhitmorel 
et al. ( |2Q11| ) also suggested this as a possible explana- 
tion lor isolated HII regions ionized by a single, massive 
star, i.e., "Single-Star" HII (SSHII) regions (see Figure 
3 from Whit more 2011). Either explanation would be 
interesting, and should be investigated in a future study. 
Figure [9] shows a similar "field" Region # 9 with primar- 
ily old red giant stars (barely discernible) , but which also 
includes a small population of isolated young blue stars 
and four compact clusters ranging in age from interme- 
diate to old. 

The right panel of Figure [71 compares the AgecMDman 
and AgecMDauto estimates for the different regions. A 
unity line is included for comparison, and shows that 
there is good agreement (7.6 a correlation) between the 
two methods which varies from the unity line by only 1 a. 
We use the objective luminosity- weighted age estimates 
(AgecMDauto) as our preferred method for estimating 
ages of the resolved stellar populations in the remainder 
of this paper. 

3.2. Comparison With Previous Age Estimates 

|Larsen et al. | ( |20TT| ) recently estimated the ages of two 
star clusters in 1V183 by comparing photometry of individ- 
ual stars in the outskirts of the clusters with isochrones. 
These clusters are within Regions #3 and #11 here. 



Larsen et al. | ( |2Q11[ ) used a more sophisticated method, 
comparing observed and predicted magnitude distribu- 
tions for red and blue stars using Hess diagrams. They 
did not however, correct the magnitudes or colors of their 
stars for the effects of spatial variation of extinction, but 
since these are intermediate age clusters with minimal 
internal extinction, this is not a large effect in these par- 
ticular cases. They find log r of 7.44 (27.5 Myr) and 7.45 
(28.2 Myr) for the clusters in Regions # 3 and # 11. We 
find similar ages of 21.9 Myr (Region #3) and 26.7 Myr 
(Region #11) for the dominant luminous stellar popu- 
lations, despite the fact that our photometry includes a 
significantly larger area around each cluster. Hence we 
fi nd good agreement be tween our age estimates and those 



oflLarsen et all (12011) 



Llach method has its own strengths and weaknesses. 
Our approach has the advantage of correcting the pho- 
tometry of individual stars for the effects of extinction to 
improve our age estimates, but the disadvantage that this 
requires imaging in more than three optical ban ds with 
a crucial need for the /7-band observation. The Larsen 



et al. (2011) method does not correct for the effects of 
extinction, but allows for a better determination of the 
star formation history over a wider age range (i.e., out 
to ~1 Gyr). 

3.3. Comparison between AgecMDauto ctnd Compact 
Cluster Ages 

Much of our recent work has involved age- dating star 
clusters using the integrated light (UBVI and Ra emis- 
sion) from the cluster and model spectral energy distribu- 
tion (SED) fitting. To compare the age of star clusters in 
each region to the age determined from individual stars 
in each region, we adopt the ages and luminositi es of 
compact clusters identified in our previous work ( |Chan-| 



dar et al. 2010). Details of the cluster age-dating method 
can be found in that paper. 

We note that the regions sampled in this study are 
large (i.e., several hundred pc 2 ) and may include stars 
and star clusters spanning a wide range of ages. This 
requires making the comparison between the ages of field 
stars and star clusters with caution. Only the largest 
star forming complexes within spiral galaxies have such 
dimensions (e.g., the giant HII region NGC604 in M33; 
Freedman et al. 2001). More typical regions are much 
smaller, and hence would not dominate the entire field. 
This will tend to weaken our correlations, especially when 
comparing resolved stellar ages with cluster ages. 

Using the ages of star clusters determined from the 
SED fitting, we compute the luminosity-weighted mean 
age of each region, as we did for stars in the same region. 
The ages of 50 regions are listed in column Ageci in Ta- 
ble [l] The upper panel of Figure [10] shows a comparison 
between the luminosity- weighted mean age (Ageci) of 
the compact clusters in a given region and our stellar 
age (AgecMDauto) estimates based on the CMDs of all 
the resolved stars in the region (as discussed in §3.1). 
The correlation between these two ages is fair (5.4 a cor- 
relation). We note that while the midpoints for the two 
methods are similar, the slope is steeper, as discussed in 
more details in the next section. 

3.4. Comparison between AgecMDauto and Age 
Estimates from Integrated Light from the Entire 
Regions 

We obtain an independent estimate of the age of each 
region by measuring the colors of the entire region, and 
perform a simple SED fit, in the same way as we did for 
star clusters (i.e., |Chandar et al. |20 10). The lower panel 
of Figure [To] shows a comparison between our age esti- 
mates for the resolved stellar populations using CMDs 
(AgecMDauto) and from the integrated light (Age Reg ) in 
each region. 

While there is a fair correlation (4.6 a) between the 
CMD and integrated photometric age estimates, there is 
also a fair amount of scatter and an apparent gap in the 
region ages in the range 6.8 < log (age/yr) < 7.2. This 
gap is similar to the well known artifact for cluster age 
estimates which is due to the looping of predicted cluster 
colors in this age range (see Chandar et al. 2010). We 
also note that the slope of the relation is steeper than 
the unity vector, with integrated light ages ranging to 
much lower ages than the CMD ages. This is similar to 
the comparison with Agecu as discussed above. This is 
probably caused by a variety of effects including: 1) the 
integrated age estimates take into account Ka emission 
(which is very sensitive to massive young stars) while 
the CMD estimates do not; 2) the 1 Myr isochrones are 
essentially on top of the 3 Myr isochrones, making it rare 
that the minimum age ever gets selected by the software 
that does the matching with the isochrones; and 3) the 
adoption of a My = — 5.5 mag cutoff removes all stars 
with ages >100 Myr from the AgecMDauto estimates, 
while the light from older stars is still included in the 
integrated light used for the AgeR eg estimate (hence, the 
AgecMDauto estimate can never be >50 Myr, while the 
AgeR e g can be older). 



This difference in slope in Figure 10 highlights the fact 
that each age-dating method has its own idiosyncrasies. 
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By comparing a number of different methods, we begin 
to understand these artifacts better, and learn how large 
the true systematic uncertainties can be. 

4. COMPARISON BETWEEN AGE C mdauto AND 
OTHER PARAMETERS THAT CORRELATE 
WITH AGE 

We are now ready to compare our age estimates for the 
resolved stellar components in 50 regions to other observ- 
ables that correlate with age. These include: (1) nu mber 



ratio of red-to-blue stas (e.g., Larsen et al. ||2011); (2) 



morphological categories (Whitmore et al. 2011) ; and (3) 
stella r surface brightness fluctuations (Whitm ore et al. | 



4.1. Comparison between AgecMDauto and Number 
Ratio of Red-to-Blue Stars 

One property that is expected to correlate with age is 
the number ratio of red-to-blue stars. At very young age, 
all stars are blue. As the population ages, the number of 
red giant stars (H-shell burning and H+He-shell burning 
stars) increases. This effect can be seen by noting that 
nearly all of the stars in Region # 12 (Figure [5| are blue, 
while there are large numbers of red stars in Regions 29 
and #2 (Figure [M. 

The number ratio of red-to-blue stars (column Red-to- 
Blue Ratio in Table [T]) is calculated by using a criteria 
that (V—/) be redder than 0.8 mag for red stars. The top 
panel of Figure pT] shows a fairly good correlation (7.2 a) 
between AgecMDauto and the Red-to-Blue Ratio. 

4.2. Comparison between AgecMDauto and 
Morphological Category 

The middle panel of Figure [TT] shows the correlation 
between our automatic CMD age estimates and the mor- 
phological classification for a given region, as defined in 
| Whit more et al. ] ( |2011| ) . Briefly, regions with Ha emis- 
sion superposed on top of the stellar component are type 
3 (emerging star clusters), regions with small Ha bub- 
bles are category 4a (very young), regions with large 
Ha bubbles are category 4b (young) , and regions with 
no Ha bubbles are category 5 (intermediate age). Based 
on this criteria, the regions displayed in Figures [5] and [6] 
are classified as categories 4a (Region # 12), 4b (Region 
#7), 5a (Region #29), and 5a/5b (Region #2), respec- 
tively. Each region was classified independently by two 
authors (HK and BCW). The mean of the two determi- 
nations is used in what follows, and is listed in column 
"Ha Morphology" of Table [Tj 

The middle panel of Figure[TT] shows that stellar CMD 
age estimates (AgecMDauto) and morphological cate- 
gories (Ha morphology) show a fair correlation (5.8 a), 
but not as good as the strong correlations (9 a for 
log r < 7 and 5 a for log r > 7) found wit h ages for 



compact clusters in Whitmore et al. ( |2011[ ). This is 



probably due to the tact that most of the regions are 
not dominated by a single age stellar population, mak- 
ing the morphological classification for an entire region 
problematic (i.e., there is strong Ha emission in part of 
the region but none in other parts). 

4.3. Comparison between AgecMDauto and Stellar 
Surface Brightness Fluctuations 



In Whitmore et al. (2011), we developed a method 
to age-date star clusters based on an observed relation 
between pixel-to-pixel flux variations (RMS) within star 
clusters and their ages. This method relies on the fact 
that the young clusters with bright stars have higher 
pixel-to-pixel flux variations, while the old clusters with 
smoother appearance have small variations. This is be- 
cause the brightest stars in a 100 Myr population have 
Mj « —3 mag, which is below our 50% photometric 
completeness limit (Mj = —3.6 mag) and well below our 
magnitude cutoff of Mj = —5.5 mag used in our stellar 
age-dating method described in §3.1. 

The bottom panel of Figure [TT] shows the correla- 
tion between the resolved stellar ages (AgecMDauto) and 
the surface brightness fluctuations (Pixel-to-Pixel RMS) 
measured in the 50 selected regions. We find little or 
no correlation (1.1 a) with a relatively large amount of 
scatter, especially for the younger clusters. Part of this 
scatter may be cau sed by the two- valued nature of the 
effect as discussed in Whitmore et al. (|2011 ), with a max- 
imum value of the surface brightness fluctuations around 
10 Myr, and lower values for both smaller and larger ages. 
Another reason for the scatter is the fact that many of 
the regions have mixed populations, rather than single 
age populations. 

5. INSIGHTS INTO THE STAR FORMATION 
HISTORY IN M83 BASED ON SPATIAL 
VARIATIONS FOR 50 REGIONS 

We first examine how the ages for the 50 regions are 
distributed in Figure [T] The boxes in this figure are color- 
coded as follows: blue (very young with ages less than 
about 10 Myr), yellow (young with ages between about 
10 and 20 Myr), and red (intermediate aged with ages 
greater than about 20 Myr). As expected, we find that 
the very young regions are primarily associated with the 
spiral arms, the young ryoungs tend to be "downstream" 
from the spiral structure (i.e., on the clockwise side away 
from the dust lanes), while the intermediate aged regions 
are found on both sides of the spiral arms in the interarm 
regions. 

A similar, more detailed analysis is being done by 
Chandar and Chien using the cluster ages in M51 (2012, 
in preparation) and M83 (2012, in preparation). They 
compare the age gradients they f ind with models devel- 
oped by |Dobbs fc Pringle] ( |2010[ ) , assuming a variety of 
different triggering mechanisms (i.e., spiral arms, bars, 
stochasticity, and tidal disturbances). 

We can perform a similar experiment here, since we 
have age estimates for each individual star. The three 
panels in Figure [12] show the distribution of stars with 
ages in the ranges 1-10, 15-35, and 40-100 Myr, respec- 
tively. The top panel is the youngest group of stars, 
clearly showing that the stars in these regions are mostly 
distributed along the active star- forming region (i.e., as- 
sociated with strong Ha emission) in the spiral arms. 
The stars in the middle panel tend to be found slightly 
downstream of the spiral arms, while the older stars are 
still farther out in the inter-arm regions, as expected. 

These diagrams can also be used to examine whether 
there is evidence that most of the young stars form in 
clust ers and clustered reg ions and then dissolve into the 
field. |Pellerin et al. | ( |2007| ) made a similar set of diagrams 
for the galaxy NGC 1313, which supported this interpre- 
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tat ion. We find that the young star samples show strong 
clustering, while the older star samples are progressively 
more uniform. Hence, these distribution maps of stars 
in the different age group support the idea that the stars 
form in clusters in spiral arms, and then mos t of the 
clusters d issolve, populating the field with stars (iLada &| 



Lada 



2003). Chandar et al. (2012, in preparation) will 



examine this subject more quantitatively in the future. 

6. SUMMARY 

Color-magnitude diagrams and color-color diagrams of 
resolved stars from the multi-band H ST /WFC3 ERS ob- 
servations of M83 have been used to measure the ages of 
stellar populations in 50 regions of this well known face- 
on spiral galaxy. The diagrams show the presence of 
multiple stellar features, including recently formed MS, 
He-burning blue- loop stars, and shell-burning red giant 
stars. Comparisons between our stellar age estimates, 
and a wide variety of other age estimators, allow us to 
investigate a number of interesting topics. The primary 
new results from this study are as follows: 

(1) An innovative new technique using a combination 
of CMD and color-color diagrams has been developed in 
order to correct for the extinction towards each individ- 
ual star and to age-date the stellar population. The mean 
extinction values for the 50 regions studied in this galaxy 
are 0.696 (A F336 w), 0.559 (A F438W ), 0.441 (A F555W ), 
0.256 (Afsiaw) mag, respectively. 

(2) The various age estimators (stellar, integrated 
light, clusters within the region) show fair correlations 
(i.e., between 4.6 a and 5.4 a). This is expected, since 
many of the regions have a mixture of different- aged pop- 
ulations within them, and each technique uses light from 
a different subset of the stars (e.g., the stellar age es- 
timates use only bright stars, while integrated light in- 
cludes all the light). A comparis on with previous age 
estimates by iLarsen et al. | (2011) shows a good agree- 
ment (see §3.2 lor details). 

(3) Comparisons between the stellar ages and other 
parameters that are known to correlate with age show a 
range from very good correlations (e.g., with red-to-blue 
ratios: i.e., 7.2 a) and fair correlations (e.g., Ha morphol- 



ogy with 5. 8 a) to little or no corre lation (e.g., pixel-pixel 
RMS from iWhitmore et al.|)2011[ ). This is expected for 
reasons similar to the ones discussed above. 

(4) The regions with ages younger than 10 Myr are 
generally located along the active star-forming regions 
in the spiral arm. The intermediate age stars tend to 
be found "downstream" (i.e., on the opposite side from 
the dust lane) of the spiral arms, as expected based on 
density wave models. A more detailed comparison with 
models with various other triggering mechanisms (e.g., 
bar s, density waves, tidal disturbance, and stochasticity; 
see |Dobbs fc Pringle |2010[ ) is in process by Chandar and 
Chien lor M51 and MS3 (2012, in preparation). 

(5) As discussed in the Appendix A, the loca tions of 



( [2005] ) 



Wolf-Rayet sources from Hadfield et al. 
broad agreement with the age estimates discussed in the 
current paper. The much better spatial resolution from 
HST shows that many of the Wolf-Rayet "stars" from 
ground-based observations are actually young star clus- 
ters. 

(6) Effects of spatial resolution on the measured 
colors, magnitudes, and the derived ages of stars in our 
M83 images are described in the Appendix B. Based on 
a numerical experiment using a star cluster NGC2108 
in the LMC, we found that while individual stars can 
occasionally show measurable differences in the colors 
and magnitudes, the age estimates for entire regions are 
only slightly affected. 

We thank Zolt Levay for making the color images used 
in Figure [T] We also thank the anonymous referee for 
helpful comments. This project is based on Early Re- 
lease Science observations made by the WFC3 Science 
Oversight Committee. We are grateful to the Director of 
the Space Telescope Science Institute for awarding Di- 
rector's Discretionary time for this program. Support 
for program # 11360 was provided by NASA through a 
grant from Space Telescope Science Institute, which is 
operated by the Association of Universities for Research 
in Astronomy, Inc., under NASA contract NAS 5-26555. 
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APPENDIX 

APPENDIX A. USING WOLF-RAYET STARS TO CHECK FOR CONSISTENCY WITH OUR AGE 

ESTIMATES 

In this paper we compared the ages resulting from a variety of age-dating methods for different star-forming regions 
in M 83. Here, we ext end our analysis by considering previously identified Wolf-Rayet stars in M83, which are ~l-4 Myr 



ol d (|Crow ther 2007), as another independent age estimate. 
Hadfield et al. I tbi)b\ ic 



" pt)05D identified 283 Wolf-Rayet rom narrow-band imaging centered on the He II 

A4686 emission line, plus spectroscopic follow-up. Seventeen of these Wolf-Rayet sources fall in our WFC3 field-of-view. 
Figure [13] shows their locations in the HST image, and Table 2 lists the coordinates of these sources, as well as the 
region in which they are located. 



The lower resolution of the ground-based Hadfield et al. (2005) study makes it difficult (in many cases) to uniquely 
identify the Wolf-Rayet sources in the HS'l image, although we can easily determine whether such a source is located 
within one of our 50 regions. Figure 13 shows a color image of M83 (center panel) with five regions, labeled A through 
E, that contain all 17 Wolf-Rayet sources. Here, the Ha and F814W images are shown in red, F555W in green, and 
F336W in blue. Black and white images are the image cutouts in the F814W and the Ha filte rs of regions A-E, as 
outlined in green in the central panel. The original finding chart of the region A taken from Hadfield et al. (2005) 
is shown in the top right panel in Figure 13 In the Ha image, the stellar continuum has been subtracted from the 
Ha image, leaving just the ionized gas. We matched the coordinates of the Wolf-Rayet sources and the HST images, 
by assuming that the Wolf-Rayet source # 99 is perfectly centered on a relatively isolated, compact Ha knot. We also 
compared the locations of the 17 matched Wolf-Rayet sources in the M83 HST images with the finding charts from 
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Hadfield et al. (2005) to verify our source matching. 

The excellent spatial resolution of the HST image also makes it possible to check whether any of the Wolf-Rayet 
sources are actually compact star clusters. While the majority do appear to be individual stars, or a dominant star 
with a little "fuzz" which is likely a faint companion star, there are at least three cases where the counterpart appears 
to be a bright compact cluster (i.e., # 74, 86, 105: circled objects in Figure 14). These three objects ar e brighter than 
My = —9 mag with the concentration index (CI) larger than 2.3, consistent with the cluster CI used in |Chandar et al.] 
(2010). There are several more cases where the counterpart is in a looser association of stars (e.g., # 78, 102, 103). 

All the Wolf-Rayet sources, except the sources in cutout E and # 108 in cutout D, are located in the spiral arms. 
This is expected, since these are the sites of most recent star formation. Nine of the 17 are in regions of strong 
Ha emission, seven more are in regions of faint Ha emission, and only one source (# 107 in cutout D) appears to have 
no Ha associated with it. 

Our primary question is whether the Wolf-Rayet sources tend to be found in regions that we estimated to have 
young ages. In Table 2 we match the Wolf-Rayet sources with the region where they are found in and find that their 
ages range from 1.1 Myr (Region #35) to 50.1 Myr (Region #41). This shows t hat such a corre l ation does exist, with 
6/8 (75 %) of Wolf-Rayet sources in regions having AgeR eg <4:.6 Myr. Similarly, |Hadfield et al. J 2005 ) found that five 
of their Wolf-Rayet sources are included in a compact star cluster catalog developed by Larseh'"(2O04), and three of 
these five clusters (#193, 179, 179) have estimated ages in the range 1.5-6 Myr, consistent with ages expected for 
Wolf-Rayet stars. 

Figure 14 shows the CMD and color-color diag ram of the measured colors from our HST observations of the object 
that best matches the Wolf-Rayet sources from Hadfield et al. (2005). We find that most of these sources are in the 
top left of the two-color diagram. We conclude that there is a strong tendency for Wolf-Rayet sources to be associated 
with the youngest regions of star formation. The good correlation provides general support for the accuracy of our 
CMD age estimates. 



APPENDIX B. EFFECT OF SPATIAL RESOLUTION 

Spatial resolution affects virtually all studies of individual stars in nearby galaxies. As is typical, in this paper we 
have assumed that the point sources identified in Section 2 are individual stars, although we also discovered a small 
fraction of pairs of nearly aligned stars, one red and one blue, that can be identified from their anomalous colors (e.g., 
"pseudo-match" Figures 2 and 4). Here, we address more generally how spatial resolution affects the measured colors 
and magnitudes, and hence the derived ages, of stars in our M83 regions. 

We begin by answering a simple question: "What would the Orion Nebula look like at the distance of M83? Would 
we be able to distinguish the four central stars that make up the Trapezium or would they appear as a single star?" 
The separation between the four stars in the Trapezium (10" at D w 0.5 kpc) would be roughly 0.001" at the distance 
of M83, approximately a factor of 40 smaller than a single pixel in our image, and hence these would appear as one 
point-like source. However, we note that the vast majority of the remaining bright stars in the Orion nebula are much 
more widely separated, and would not suffer from this problem. We also note that our selection criteria are designed 
to avoid including stars in the crowded central regions of star clusters and associations, which also minimizes the 
problem. 

We perform a more quantitative analysis by degrading HST /AC S images of the compac t star cluster NGC2108 in 
the Large Magellanic Cloud (LMC) (ID:10595; PI: Goudfrooij, see Goudfrooij et al. 2011), which is at a distance of 
~ 50 kpc, by a factor of 100 (i.e., equivalent with the distance of ^ 5Mpc), as roughly appropriate for M83. The 
degraded images, created from a combination of rebinning and Gaussian smoothing to approximately mimic the point 
spread function of ACS, are then run through our normal object-finding software described in Section 2. The numbers 
of objects found in the degraded NGC 2108 images are ~20 in the F555W and F814W ACS images. Figure 15 shows 
the F555W ACS image of NGC 2108, along with the degraded images and a color image (i.e., a color composite of 
F435W, F555W, and F814W ACS images). 

Additionally, we then perform aperture photometry on the detected source positions in the non-degraded ACS 
images using a radius 300 pixels to collect light from all stars that would fall within our aperture if NGC 2108 was 
located at a distance of 5Mpc (i.e., M83). Objects along the edge and within 1000 pixels of the center of the cluster in 
NGC 2108 are discarded, as would be the case for the corresponding photometry in M83. Figure 16 shows the extent 
of the 300 pixel aperture for several sources in yellow, as well as the nearest dominant star in red. The magnitude 
of the dominant stars (using a 3 pixel aperture and the appropriate aperture correction) is considered the "truth" 
measurement while the magnitude in the corresponding 300 pixel aperture is the value that would be determined at 
the distance of M83. We calculate the magnitude and color offsets, i.e., Al and A(V — I) mag of the dominant stars, 
by taking the difference between the 3 pixel and 300 pixel aperture magnitudes. 

One complication is that our default sky subtraction method would remove the vast majority of faint and moderate 
brightness stars in the background annulus in the NGC 2108 image, but these stars would not be detected individually 
at the distance of M83 and hence would be included in the sky measurement. We therefore adjusted the "clipping" 
parameters in the sky subtraction algorithm to remove only the bright stars in NGC 2108, thereby mimicking the case 
for M83 to the degree possible. 

In Figure 17 we show the original CMD for Region 3 in M83 (solid points), compared with a version where the 
magnitudes and colors have been 'perturbed' by an amount Al and A(V — I) determined by matching to the closest 
values of (V—I) from the NGC 2108 data experiment (open points). As expected, corrections for the degraded spatial 
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resolution in M83 tend to make the corrected magnitudes slightly fainter and redder, as can be seen in Figure 17 
The mean difference in V is —0.24 ± 0.30 mag with a range from +0.20 to —0.65 mag. In /, the mean difference is 
—0.40 ± 0.39 mag with a range from +0.07 to —1.18 mag. We note that the largest difference is seen for object 10, as 
shown in Figure 16 , where several relatively bright red stars in the 300 pixel aperture cause the largest correction (i.e., 
— 1.18 mag in I and redward by 0.63 mag in (V — /)). Object 40, dominated by a single very bright star, shows the 
opposite extreme, with a change of —0.02 mag in V and —0.01 mag in I. We also note that object 31 does not make it 
into our sample since the two bright stars are far enough apart that it is clear from the degraded image in Figure 15 
that this is not a single star. This object gets removed in the DoPHOT photometry as described in § 2, and hence is not 
used in the experiment described above. 

We then rerun our age-dating software on the corrected CMD in Figure [l7| We find a mean age of 23.8 Myr compared 
to the original value of 21.9 Myr. Hence, while a specific "star" may be affected by a sizeable amount, the overall 
effects of degraded spatial resolution are relatively minor for our study. 

To summarize, we have performed a numerical experiment using observations of an intermediate-age star cluster in 
the LMC (NGC 2108) for our "truth" image. Using a factor of 100 in spatial degradation, roughly appropriate for 
our measurements in M83, we find that the typical corrections to our photometry are on the order of a few tenths 
of a magnitude, although much larger values are possible in specific cases. We conclude that while spatial resolution 
can result in measurable differences in the luminosities and colors of bright stars at the distance of M83, it does not 
strongly affect the age estimates for the 50 regions studied in this paper. 
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Fig. 1. A color composite of the M83 WFC3 images [Image Credit: Zolt Levay (STScI)]. The F336W image is shown in blue, the F502N 
([O III]) image in cyan, the F555W image in green, and the combined F814W and F657N (Ha) image in red. The 50 selected regions are 
outlined in boxes of three different colors based on the values of AgecMDauto determined in this paper (see §3.1 for details). The regions 
with ages in the range 1-10 Myr are outlined in blue, ages of 10-20 Myr in yellow, and ages greater than 20 Myr in red. 
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Fig. 2. — Image cut-outs (15'/2xl3'/9, 340pcx310pc) of Regions # 26 and # 27 in the F555W (top) and the color-composite images 
(bottom). Stars circled in red have (V — I) >0.8 mag. The "postage stamp" images on the left are"pseudo" matched candidates, while the 
images on the right appear to be heavily reddened (see §2.2 for details). 
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Fig. 3. — Top: The color-magnitude diagram (CMD) and color-color diagram of all stars in our M83 image corrected for Galactic 
foreground extinction but not for internal extinction. Padova isochrones of logr (age/yr) = 6.55, 6.85, 7.0, 7.2, 7.5, 7.8, and 8.0 for a 
metallicity of Z=0.03 (1.5 Zq) are overlaid in both panels. The dashed line in cyan represents the 50% photometric completeness level. 
The arrow in each panel indicates the direction of the reddening vector with Ry =3.1. Magnitudes are on the Vega scale. The distance 
modulus (m — M)o = 28.32 mag is used to calculate the absolute magnitude Mp8i4W- Bottom: Histograms of the distribution of ages and 
masses, as determined from the CMD (see §3.1 for details). The age quoted in the bottom left panel is the luminosity- weighted mean age. 
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Fig. 4. — Top; The CMD and color-color diagram of all stars corrected for the internal extinction determined for each individual star. The 
blue and two red dotted lines indicate boundaries for the locations of stars in the color-color diagram uncorrected for internal extinction (see 
§2.4 for discussion). Two arrows in the color-color diagram show reddening vectors: one with the standard reddening vector (Ry — 3.1) 
and the other with the flatter reddening vector (Ry = 5.7) which appears to be more appropriate for some of the data in the green triangle 
(see §2.4 for discussion). The dashed line in magenta shows the Mj<— 5.5 magnitude cutoff used for the AgecMDauto estimate. Bottom: 
Same as described in Fig. [3] 
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Fig. 7. — Correlation plots between the ages of the stellar populations in the 50 regions from CMD isochrone fitting. (Left): Correlation 
between the manual ages determined by HK and BCW. (Right): Correlation between the average of manual ages and automatic ages (see 
§3.1 for details). The best linear fits are shown in black with the significance of the correlation in unit of a (i.e., slope/uncertainty) in the 
top left of each panel. The dashed line is the unity line. 
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Fig. 8. — Histograms of the distribution of ages and masses of the individual stars in Regions # 12, # 7, # 29, and # 2, as determined by 
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Fig. 9. — Image cutout (18 f /Sx 16'/9, 410pcx380pc) of the area near Region 9 as an example of a "field" region with old red giant stars, 
and a small population of isolated young blue stars. See §3.1 for for detail. 
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Fig. 10. — (Top): Correlations between the ages of stars (AgecMDauto) and clusters (Ageci)- (Bottom): Correlation between the ages 
of stars determined by AgecMDauto and the integrated light within the region (AgeR eg ). See §3.4 for details. The best linear fits are 
shown in black solid line with 5.4 a (top) and 4.6 a (bottom) correlations. The dashed line is the unity line. 
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Fig. 12. — Spatial distribution of stars with ages of 1-10, 15-35, and 40-100 Myr determined in this study. Each panel has the same 
orientation and scale as Figure 1. The top panel is the youngest group of stars, clearly showing that the stars in these regions are mostly 
distributed along the active star-forming region (i.e., associated with the strong Ha emission) in the spiral arms. The stars in the middle 
panel tend to be found slightly downstream of the spiral arms while the older stars are still farther out in the inter-arm regions, as expected. 
See §5 for details. 
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Fig. 13. — Distribution of Wolf-Rayet sources in our HST WFC3 M83 images. The color image at the center is the UBVIHa composite 
and the image cut-outs are the Ha {left or top) and F814W {right or bottom) images for regions A-E. The top right panel is an original 
finding chart (Hadfield et al. 2005 A4684 narrow-band FORS2 image) for the region A. Numbers in blue are the IDs of Wolf-Rayet sources 
identified by |Hadfield et al. Hi2U05) . 
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Fig. 14. — The CMD and color-color diagram of Wolf-Rayet sources in M83 marked in Figure fl3] Three sources (#78, 86, 105) circled 
in both diagrams are identified as clusters in our M83 WFC3 images. A twice solar metallicity BC09 cluster model (Bruz ual Sz Chariot 
2009, private communication; see also Bruzual & Chariot 2003) is shown in the right panel (dots), and Padova stellar models ( |Marigo et al. | 
2008) are shown in both panels. No correction has been made for internal reddening, which appears to be minimal for all but one object. 
The arrow indicates the Galactic reddening vector. 




Fig. 15. H ST/ ACS images (290" x290", 68pcx68pc) of the compact star cluster NGC2108 in the LMC: a color composite of the 
F435W, F555W, and F814W ACS images (left), the F555W ACS image (middle), and the degraded F555W ACS image (right) by a factor 
of 100 (from a distance of ~50kpc to a distance of ~5Mpc). See Appendix B for details. 
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Fig. 16. — Image cutouts of H ST/ ACS F555W image of NGC 2108. The large circles in yellow show the 300-pixel aperture (corresponding 
to a ~3 pixel aperture at the distance of M83) centered at the location determined from the degraded image in Figure 15. The location of 
the dominant star (used to determine the "truth" value) is shown in red. See Appendix B for details. 
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Fig. 17. — The original CMD of Region 3 in M83 (solid points) with the effect of simulation representing a factor of 100 degration in 
spatial resolution (open points). The AgecMDauto of this region changes from 21.9 Myr (determined from the solid points) to 23.8 Myr 
(from the open points). The difference in the age estimates is within ~ 1 a error listed in Table 1. See Appendix B for details. 



RESOLVED STARS IN M 83 



27 



cocNLO^cNi^o^b-ococN^OLOOCN^coLOCNb-^ocoi^b^ 



6 

CD 

6 
13 



bO 

CD ^ 



O 

o c 

CN 



31 

o 

"Sb 



CD 


T— 1 




b- 




CN 




oo 


1—1 






co 


OS 


co 


co 

OS 


oo 


o 


CD 


oo 


b- 


CD 


oo 


OS 


oo 


LO 

co 




b- 


CD 


co 


co 




oo 


LO 


LO 


b- 


OJ 


b- 


b- 


O 




co 


CD 


b- 




o 


CN 

oo 


o 


O 


co 


T— 1 


CO 


b- 


co 


CD 


OS 


CN 




b- 


co 


CN 


o 


oo 


b- 


oo 


C5 


CN 




o 


CN 


i—l 


CN 


LO 


LO 


O 




b- 


CN 




O 


oo 


CN 




CN 


oo 


OS 


C5 


CN 


co 


CN 


CM 


CN 


co 


T— 1 




CN 


T— 1 


i—l 


b- 


i—i 


i—i 


b- 


co 






1—1 


CN 


co 


CN 






CM 




i—l 


i—i 


CN 




i—i 


co 




i—l 


co 


CN 


i—l 








1—1 


CN 


i—i 





oo 




i—i 


oo 


b- 


OS 


LO 


OS 


CN 


CD 


CD 


oo 


oo 


LO 


CN 


oo 


OS 




CN 


CN 


co 


i—i 


o 


CD 


b- CM 


b- 


b- 


co 


b- 


b- 


CD 


b- 


co 


CN 


CD 


CN 


o 




CN 


OS 


CN 


o 


LO 


OS 


CD 


b- 


CD 


b- CM 


oo 


1—1 


CD 


LO 


b- 


oo 


o 


CD 




CN 


oo 


LO 


b- 


b- 


CN 


o 


co 


o 


OS 


CD 


o 


o 


i—i 


CN 


CM Ttf 


b- 


i—i 




i—i 


oo 


o 


i—i 


OS 


oo 


oo 




CN 


b- 


b- 




co 


o 




OS 




CN 


o 


O LO 


CN 


CN 




oo 


CD 


CN 


co 


co 


LO 


b- 


1—1 


b- 


co 




co 






CN 


co 


00 


co 


oo 


LO 


oo 


CD CO 


LO 


q 


CN 


co 




co 


CN 


co 


b- 


CN 


CN 


LO 


i—i 






CN 


co 


CN 


oo 


oo 






LO CD 


O 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d d 


d 


i — 5 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d d 



LO 

co 


LO 

co 


o 
q 


o 
o 


o 
o 


LO 
CN 


o 

LO 


LO 

co 


o 

LO 


o 

CD 


o 
q 


o 
i—i 


o 

LO 


LO 
CN 


o 
q 


o 
i—i 


LO 
CN 


LO 
OS 


LO 

CD 


O 
O 


O 

q 


o 

CD 


LO 
OS 


o 
o 


o o 

CD LO 


o 

LO 


lO 
CN 


o 
q 


lO 

co 


o 

lO 


o 
q 


o 
i—i 


lO 

CD 


O 
i—l 


o 

1—1 


lO 

co 


o 
i—i 


o 
q 


o 

lO 


o 

oo 


o 
q 


o 

CN 


lO 
OS 


lO 

b- 


o 

1—1 


o 
q 


o 
q 


CM CO 


LO 


LO 


LO 




LO 


LO 






LO 




LO 




LO 


LO 


LO 




LO 








LO 








CO TJH 






lO 


lO 




lO 


lO 






id 


id 




id 






id 


id 








id 


id 


id tjh 



co 


CN 


b- 




b- 


o 


CD 


o 


CN 


CD 


i—i 


CD 


CN 




lO 


o 


b- 


CD 


CN 


oo 


b- 


lO 


o 


OS 


o co 


b- 




lO 


co 




oo 


CN 


oo 


b- 


b- 


lO 


o 


co 


OS 


co 


lO 


lO 


OS 


OS 


o 


co 


CN 


CO CD 


OS 


OS 


i—i 


b- 


CD 


o 


o 


o 




lO 


i—i 


co 


co 


CD 


co 


o 


b- 




CD 


oo 


C5 


CD 


oo 


b- 


O CD 


lO 


i—i 


CN 


oo 


o 


i—i 


i—l 


i—i 


b- 




CD 


i—i 


co 


CN 


lO 


CN 


oo 


CD 


b- 


o 


co 




CO OS 




b- 


oo 


i—i 


CN 




co 


CN 




1—1 


OS 


1—1 


CD 


co 


i—i 


CN 


co 




CN 


i—i 


CN 


1—1 


o 


o 


o ^ 


1—1 


CN 


CD 


lO 




CD 


oo 




co 


OS 


oo 


CD 


b- 






oo 






i—i 


CD 


co 


CN 


^ CM 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 




d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d d 



CN 


OS 


OS 


CN 


CD 


lO 




oo 


OS 


oo 


OS 




co 


CN 


CD 




CN 


oo 


o 


CN 


OS 


oo 


OS 


o 


b- i-H 


CN 


CN 


OS 


CN 


CN 


CN 


o 


o 


CD 


o 


co 


b- 


lO 




CN 


o 


CD 




CN 




lO 


OS 


CO OS 


oo 


co 


b- 




CD 


oo 


OS 


i—i 




co 


co 


oo 


oo 


CD 


oo 


co 


lO 


CN 


o 


oo 


co 




b- 


OS 


05 OS 


co 


i—l 


lO 


co 


co 


co 


OS 


lO 


b- 


lO 




CN 


CN 


oo 


lO 


OS 


1—1 


b- 


i—l 


o 


lO 


i—i 


CO CD 


lO 


lO 


CD 




io 


lO 


oo 


b- 


co 


CD 


CD 


CD 


co 


lO 


CD 


lO 


b- 


05 


05 


CD 


OS 


oo 


C5 


CD 


TJH OS 


o 


CD 


CD 


CD 


CD 




CD 


b- 


lO 






b- 


oo 


CD 


b- 


CD 


lO 


lO 


lO 








CO TJH 


q 


o 


O 


o 


O 


o 


o 


o 


o 


o 


o 


o 


o 


O 


O 


o 


o 


O 


O 


O 


o 


o 


o 


o 


O O 




O 


O 


o 


o 


o 


o 


o 


o 


o 


o 


O 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o o 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d d 



co 


b- 


CN 


oo 


OS 


OS 


co 


b- 


CD 


CN 


2 


q 


OS 


CN 


b- 


co 








CD 


co 


co 


^N 


LO i-H 


b- 


b- ^ CN 


co 


co 




O 


OS 






i—i 


oo 


CN 




o 


CD 


o 


o 


CN 






CN 


CD 








d 


i-H 


co 


CO 


co 


lO 


id 


CN lO 


CN 


^ CN LO 




i—i 




CN 






b- 


lO 


CN 


CN 


co 




CN 




i—l 




co 




CN 








CN 


CN 




in io 




^ b- CN 




LO 


CD 


i—l 







ONCDO 

TJH LO TJH TJH 



^ O ^ 9 

CO 4 o CN b- 



OS 




b- 


9 


O 


b^ 


b^ 


oo 


CD 






i—l 


CN 


CN 



^CN^^^OO^OOCD 



CN CD ^ 
CO ^ ^ 



COCNb-LOOi^CNLOO^CDOCOCOO^LOCDOCOLO^LOCDOOLOOi^^CDi^OOCNCO 

i-5cNi-5di-5i-5i-5dcNCNi^di-5dddcNddddddd^ 



^Loqoq^qq^^i 



.co^hocooonoo^hoolooo^. 



,0100010 10, 



: O LO CO CN m i-H , 



o 


LO 


O 


o 


LO 


o 


LO 


LO 


o 


o 


LO 


o 


o 


LO 


o 


o 


LO 


o 


LO 


o 


CN 
i-H 


b- 

CN 


OS 
1—1 


b- 

1—1 


CN 
CN 


LO 

CN 


CN 


o 

1—1 


LO 

CN 


LO 

CN 


o 

1—1 


b- 

1—1 


o 

1—1 


LO 
1—1 


OS 

1—1 


LO 

CN 


1—1 


CD 


i—l 
i—l 


1—1 



2^ 



o 


o 


O 


o 


o 


o 


O 


o 


o 


o 


o 


o 


O 


o 


o 


o 


O 


o 


o 


o 


O 


O 


O 


O 


o o o 


o 


O 


o 


o 


o 


o 


o 


o 


o 


O 


o 


o 


o 


o 


o 


O 


o 


o 


o 


O 


o 


o o 


1—1 


OS 


i—l 


CD 


OS 


1—1 


LO 


CD 


oo 


b- 


LO 


o 


oo 


i—i 


oo 


o 


co 


C5 


oo 


oo 


co 


CN 


i—l 


LO 


Oi CD CD 


co 


o 


1—1 


i—i 


oo 


b- 


co 


co 


OS 


i—i 


i—i 


i—i 


oo 


o 


b- 


b- 


CD 




OS 


i—i 


OS 


CN CN 


co 


CN 


co 


CN 


CN 


CN 


CN 


CN 


CN 




CN 


CN 




co 


i—i 


co 


CN 


1—1 


i—i 




CN 


co 




CN 


i — 1 CN i — 1 


co 


CN 


co 


LO 


co 


CN 


co 


CN 




LO 


CN 






co 


co 


CN 


CN 




CN 


co 


CN 


CN CN 


X 


X 


X 


X 


X 


X 


X 


X 


X 


X 


X 


X 


X 


X 


X 


X 


X 


X 


X 


X 


X 


X 


X 


X 


XXX 


X 


X 


X 


X 


X 


X 


X 


X 


X 


X 


X 


X 


X 


X 


X 


X 


X 


X 


X 


X 


X 


X X 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o o o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o o 


1—1 


LO 




LO 




oo 


LO 


LO 


CN 






LO 






o 


^ 


co 


co 


co 


oo 


CN 


CN 




05 


oo oo O) 


CD 


oo 




b- 




b- 


o 


CD 


OJ 




oo 


CD 


CN 


o 


co 


LO 




OJ 


CN 


LO 


b- 


CO CN 


co 


co 


co 


CN 


CN 


CN 


CN 


CN 




CN 


CN 


CN 


co 




CN 


CN 


CN 


CN 




1—1 


co 


CN 


1—1 


i—l 


i-H i-H i-H 


1—1 


1—1 


CN 


CN 


CN 


co 


co 


CN 


co 


co 


1—1 


CN 


CN 


CN 


CN 


CN 




i—l 


CN 


CN 


CN 


CO CN 


OS 


o 


CN 


LO 


o 


CN 


CN 




oo 


oo 




oo 


oo 


CD 


^ 




CD 


oo 


oo 


b- 


b- 


b- 


CD 


o 


CN b- CD 


CD 


b- 


CN 


o 


oo 


OS 


1—1 


oo 


oo 




OS 


oo 


OS 




o 


CN 


o 




CD 


b- 


co 


b- co 


1—1 


CO 


b- 


co 


CN 


co 


b- 


CN 


oq 






oq 




OS 


o 




LO 






o 


oq 


oo 


b- 


oo 


LO CO LO 


b- 




co 


CN 


CN 


b- 


co 




q 


oq 


co 




LO 


q 


CN 


05 


oo 


CN 


LO 


i—l 


oo 


CN CN 


d 


1—1 


d 


oo 


CN 


CD 


oS 


CN 


i—i 


co 


b^ 


oo 




LO 


d 


CN 


cd 


CD 


b^ 




CD 




d 


b^ 


d oo cn 


co 






CD 


CN 




d 




co 




d 










cd 


d 


d 


CN 




d 


d d 


co 


CN 


i—i 






LO 


q 


CN 


^ 




LO 


i—i 


CN 


1—1 






q 


LO 


i—i 


CN 


CN 


CN 


co 


co 


LO ^ LO 


LO 


1—1 


i—i 


i—i 


CN 






co 




LO 


LO 


o 


o 


i—i 


i—l 


1—1 


o 






co 


CN 


th q 


CN 


CN 


CN 


CN 


CN 




CN 


CN 


CN 




i—l 


CN 


CN 


CN 


i—i 




CN 


1—1 


CN 


CN 


CN 


CN 


CN 


CN 


CN CN CN 


CN 


cd 


cd 


co 


co 


cd 


cd 


cd 


cd 


cd 


cd 














co 


co 


cd 


co 


co co 


LO 


LO 


LO 


LO 


LO 


LO 


LO 


LO 


LO 


LO 


LO 


LO 


LO 


LO 


LO 


LO 


LO 


LO 


LO 


LO 


LO 


LO 


LO 


LO 


LO LO LO 


LO 


LO 


LO 


LO 


LO 


LO 


LO 


LO 


LO 


LO 


LO 


LO 


LO 


LO 


LO 


LO 


LO 


LO 


LO 


LO 


LO 


LO LO 


d 


6S 


6S 


6S 


6S 


d 


6S 


6S 


6S 


6S 


6S 


6S 


6S 


6S 


6S 


6S 


6S 


d 


d 


d 


d 


d 


d 


d 


cji cji cji 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d d 


CN 


CN 


CN 
i 


CN 
i 


CN 


CN 
1 


CN 
i 


CN 


CN 


CN 
i 


CN 


CN 


CN 
i 


CN 
i 


CN 


CN 
i 


CN 
i 


CN 


CN 


CN 
i 


CN 
i 


CN 


CN 
i 


CN 
i 


CN CN CN 
i i i 


CN 


CN 


CN 
i 


CN 
i 


CN 


CN 
i 


CN 
i 


CN 


CN 


CN 
i 


CN 


CN 


CN 
i 


CN 
i 


CN 


CN 
i 


CN 
i 


CN 


CN 


CN 
i 


CN 


CN CN 


oo 


i-H 


O 


oo 


LO 


oo 


b- 


oo 


CN 


LO 


CN 


O 


CN 


i—l 


co 


CN 


i—l 


1—1 


i-H 


CN 


LO 


OS 


b- 


i—l 


CN 00 CO 


OS 


i-H 


oo 


b- 


CN 


b- 


LO 




oo 


co 


O 


OS 


OS 


O 


i-H 


CN 


b- 


b- 




OS 


CD 


OS ^ 


b- 


i-H 






co 


1—1 


O 


LO 


CN 


i—l 


b- 


oq 




co 


LO 


co 


co 


CN 


OS 


CD 


LO 




q 


q 


^ co b- 


o 


co 




LO 


q 


oq 


CD 


LO 


CD 


OS 


oq 


CN 


co 


CO 




OS 


oq 


b- 


q 


OS 


00 


b- CN 


C5 


i—l 


CN 


CN 


co 




^ 




LO 


LO 


LO 


LO 


CD 




b^ 


oo 


oo 


d 


oo 


oo 


d 


d 




d 


oo i-5 d 


d 


i—i 


d 


d 


oo 




d 


d 


b^ 


oo 


d 


d 


d 


oo 




CN 


CN 


i—l 


i—i 


CN 


i—i 


d t-h 


LO 


O 


O 


O 


o 


o 


o 


o 


o 


o 


o 


o 


O 


o 


o 


o 


o 


o 


o 


o 


o 


i—i 


i—i 


o 


O i-H i-H 


i—i 


i—i 


i—i 


o 


o 


i—i 


1—1 


o 


o 


o 


o 


i—i 


o 


o 


O 


O 


O 


O 


o 


O 


o 


o o 


CD 


b^ 










t"^ 


r"^ 


b^ 


t"^ 


r"^ 


b^ 


t"^ 




b^ 




















b-^ b-^ b-^ 


















r"^ 












b^ 




t"^ 


b^ 






r"^ 




CO 


co 


co 


co 


co 


co 


co 


co 


co 


co 


co 


co 


co 


co 


co 


co 


co 


co 


co 


co 


co 


co 


co 


co 


co co co 


co 


co 


co 


co 


co 


co 


co 


co 


co 


co 


co 


co 


co 


co 


co 


co 


co 


co 


co 


co 


co 


co co 


cd 


cd 


cd 


cd 


cd 


cd 


cd 


cd 


cd 


cd 


cd 


cd 


cd 


cd 


cd 


cd 


cd 


cd 


cd 


cd 


cd 


cd 


cd 


cd 


cd cd cd 


cd 


cd 


cd 


cd 


cd 


cd 


cd 


cd 


cd 


cd 


cd 


cd 


cd 


cd 


cd 


cd 


cd 


cd 


cd 


cd 


cd 


cd cd 


1—1 


i—i 


i—i 


i—i 


i—i 


i—i 


i—i 


i—i 


i—i 


i—i 


i—i 


i—i 


i—i 


i—i 


i—i 


i—i 


i—i 


i—i 


i—i 


i—i 


i—i 


i—i 


i—i 


i—i 


i-H i-H i-H 


i—i 


i—i 


i—i 


i—i 


i—i 


i—i 


i—i 


i—i 


i—i 


i—i 


i—i 


i—i 


i—i 


i—i 


i—i 


i—i 


i—i 


i—i 


i—i 


i—i 


i—i 


i-H i-H 




















o 


i—i 


CN 


co 




LO 


CD 


b- 


oo 


a 


o 




CN 


co 




LO CD b- 


oo 


a 


o 




CN 


co 




LO 


CD 


b- 


oo 


a 


o 




CN 


co 




LO 


CD 


b- 


oo 


OS O 


1—1 


CN 


co 




LO 


CD 


b- 


oo 


OS 


i—i 


i—i 


i—l 


i—i 


i—i 


i—l 


i—i 


i—i 


i—i 


i—i 


CN 


CN 


CN 


CN 


CN 


CN CN CN 


CN 


CN 


co 


co 


co 


co 


co 


co 


CO 


co 


co 


co 




















^ LO 



73 



I 



• co 

CO 

^ co 

coo 

V> v 

^ CO 



3 bb 

CO 4f 

cooca 
co H 

co • CO 

9 Q ^ 

J H rD ' 

^ bJO 0) 
?h d h 

^1 . 

^ ^ % . 

00 73 ^ , 

CD ct3 bO 
CD co CD 
CO co"'^ ^ • 

c6 ^ 13 S 

© H ,H 

*d .t| 73 co 

1 ij i- 

CD CD O 

^ 5 CD CD 

cd ^ bO > 

§ ^-2 S 
e 8 ^5 

cd u a 
| g>S* 

h H h 

CD CO co 3 ■ 

sail 
3 s l § 

CD . S 
^ _^ CD ^ 

° § CD 3. 
S ^ CO 

^ S fl O 1 

S « d § 

£ 8,-8 1' 

.2 bO > *w 
CD 7 O 



S "3 CD ' 

"i CD ■ 

^ CD CD > 

< ^ H ^ 



2 ^ 

CD bO 



Kim et al. 





Table 2. Wolf-Rayet Star 


Candidates in 


M83 






RA 


DEC 


n aaj xeia 


Reg# 


Aqe Rea (Myr) 


X 


Y 


Comments 


13:37:00.41 


-29:52:54.1 


72 






3874 


606 


faint Hen 


13:37:01.18 


-29:52:53.9 


73 






3618 


609 


faint Ha 


13:37:01.42 


-29:51:25.8 


74 


4 




3512 


2821 


strong Ha 


13:37:03.06 


-29:50:49.5 


78 


7 


4.2 


2982 


3698 


faint Ha 


13:37:04.65 


-29:52:47.5 


82 






2710 


772 


strong Ha 


13:37:04.65 


-29:50:58.4 


86 


12 


2.4 


2477 


3505 


strong Ha 


13:37:07.35 


-29:49:39.2 


95 


20 


0.8 


1513 


3306 


strong Ha 


13:37:07.54 


-29:52:53.6 


96 


41 


50.1 


1552 


620 


strong Ha 


13:37:07.96 


-29:52:07.9 


99 






1418 


1762 


faint Ha 


13:37:08.25 


-29:51:13.6 


100 


24 




1311 


3116 


faint Ha 


13:37:08.40 


-29:52:54.9 


102 







1266 


588 


strong Ha 


13:37:08.53 


-29:52:12.0 


103 


35 


1.1 


1224 


1659 


strong Ha 


13:37:08.70 


-29:52:28.9 


105 


38 


2.1 


1173 


1239 


strong Ha 


13:37:08.91 


-29:52:05.7 


106 


31 


4.4 


1101 


1818 


faint Ha 


13:37:09.02 


-29:52:45.2 


107 


39 


24.0 


1067 


825 


no Ha 


13:37:09.80 


-29:52:36.2 


108 






820 


1049 


strong Ha 


13:37:10.42 


-29:51:28.0 


109 


26 


4.6 


610 


2756 


faint Ha 



